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Abstract 

Interfaces in a model with a single, real nonconserved order pa- 
rameter and purely dissipative evolution equation are considered. We 
show that a systematic perturbative approach, called the expansion 
in width and developed for curved domain walls, can be generalized 
to the interfaces. Procedure for calculating curvature corrections is 
described. We also derive formulas for local velocity and local sur- 
face tension of the interface. As an example, evolution of spherical 
interfaces is discussed, including an estimate of critical size of small 
droplets. 
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1 Introduction 



Important aspect of dynamics of phase transitions in condensed matter is 
time evolution of an interface separating a retreating phase from the new 
one. Studied in the framework of Ginzburg-Landau type effective macro- 
scopic models the interface can be regarded as a kind of smooth, asymmetric 
domain wall subject to a transverse force. The asymmetry and the force are 
due to a difference of potential energy across the interface. Pertinent evolu- 
tion equations for order parameters typically are nonlinear partial differential 
equations. In general they imply rather nontrivial phase ordering dynamics, 
see, e.g., review article Relativistic version of the problem, not consid- 
ered here, also is interesting because of its connection with field-theoretical 
cosmology, 0. 

Recently, evolution of ordinary domain walls has been studied with the 
help of Hilbert-Chapman-Enskog method applied in a suitably chosen comov- 
ing coordinate system ||, ||, @ - a systematic and consistent perturbative 
scheme has been developed. It yields the relevant solutions of the evolution 
equations in the form of expansion in a parameter / which can be regarded as 
a measure of width of a static planar domain wall. Consecutive terms in this 
expansion contain extrinsic curvatures of a surface comoving with the wall, 
and also certain functions (below denoted by Ck) which can be regarded as 
fields defined on that surface and coupled to the extrinsic curvatures. In the 
present paper, which is a sequel to the previous paper |J, we show that that 
perturbative expansion can be generalised to the case of curved interfaces. 

The Hilbert-Chapman-Enskog method and the comoving coordinates tech- 
nique, which we have learned from [g], 0, respectively, have already been 
used in theoretical investigations of planar interfaces 0, and of curved ones 
in superconducting films @. We apply these tools to curved interfaces in the 
three dimensional space, in a Ginzburg-Landau type model defined by formu- 
las (1), (2) and Eq.(3) below. Interfaces in this model have been investigated 
for a long time, see, e.g., |J. Our main contribution consists in providing a 
systematic iterative scheme for generating the relevant solution in the form 
of a perturbative series. The role of small parameters is played by the ratios 
lo/Ri, where Iq is the width and Ri curvature radia of the interface. In spite 
of nonlinearity of the evolution equation the perturbative contributions can 
be generated in surprisingly simple manner. This is achieved by introduc- 
ing the functions Ck which saturate certain integrability conditions. As an 
application, we derive a formula for local velocity of the interface with curva- 
ture corrections included, and we discuss critical size of nucleating spherical 
droplets. 

We consider a system described by a real, scalar, non-conserved order 
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parameter with the free energy F of the form 



J \2 dx a dx a v ; 



where 

V = A$ 2 + 5$ 3 + C$ 4 . (2) 
Time evolution is governed by the dissipative nonlinear equation 

T ^M = KA*-V'(*). (3) 

Here (x a ) Q , = i 2 ,3 are Cartesian coordinates in the space, V denotes the deriva- 
tive dV/dQ, and K, 7, A, B, C are positive constants. The free energy of the 
form (1) arises in, e.g., de Gennes-Landau description of nematic-isotropic 
transition in nematic liquid crystals in a single elastic constant approxima- 
tion (L 2 = 0) ||lQf . Then $ = corresponds to the isotropic liquid phase, 
while in the nematic phase $ 7^ 0. 

The concrete form (2) of the potential has the advantage that solution 
of Eq.(3) describing a planar interface has a simple, explicitly known form. 
It can be found in, e.g., ||. The planar interface plays important role in 
the perturbative scheme: the main idea is that there exist curved interfaces 
which do not differ much from the planar one if considered in appropriately 
chosen coordinate system (which in particular should comove with the inter- 
face). Therefore, one may hope that $(5f, t) for such curved interfaces can 
be calculated perturbatively, with the planar interface giving the zeroth or- 
der term. To this end, it is necessary to introduce the comoving coordinate 
system explicitly, and to give a prescription for the iterative computation of 
the perturbative corrections. 

The plan of our paper is as follows. In Section 2 we recall the planar 
interface and the comoving coordinates. This Section contains preliminary 
material quoted here for convenience of the reader as well as in order to 
fix our notation. In Section 3 we describe the perturbative scheme for the 
curved interfaces. In Section 4 we present formulas for the velocity and the 
free energy of the curved interface. Section 5 is devoted to a discussion 
of spherical droplets of the stable phase which nucleate during the phase 
transition. Several remarks are collected in Section 6. In Appendix A we 
construct solutions of linear equations obeyed by corrections to transverse 
profile of the interface. Appendix B contains a brief discussion of stability of 
the interface. 
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2 The preliminaries 

2.1 The homogeneous planar interface 

Let us assume that the planar interface is perpendicular to the z axis [z = x 3 ) 
and homogeneous. Then $ depends only on z and t, and Eq.(3) is reduced 
to 

7 <9 t $ = Kd 2 z <& -V'(§). (4) 

The interface type solution interpolates smoothly between minima of 

V when z changes from — oo to +oo. V given by formula (2) has two minima 



where Iq is given by formula (9) below. The corresponding phases we shall 
call isotropic and ordered, respectively. Let us multiply Eq. (4) by d z Q and 
integrate over z. The resulting identity 
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/+oo 
dzd z $d t <$> = V($_) - V($+), 
-oo 



implies that d t & ^ if the minima are nondegenerate. Further assumption 
that the interface moves in a uniform manner with velocity Vq, that is that 

$(z,t) = $(Z - ZQ-V Q t), 

leads to the formula 

/+oo 
rfZ$ 2 (Z) = V($ + )-V($_), (5) 

where 

Z = z — zo — vot. (6) 

Hence the interface moves towards the region of higher potential V($±), as 
expected. 

The exact solution of Eq.(4) has the following form [|] 



$0 = -J2L I (7) 

V 8l 2 C 1 + exp(-Z/2/ ) ' K) 



where 



and 



v 



9K 



B - V9B 2 - 32AC 



V 32 7 2 C 

lo 1 = ( 35 + V9B^32AC) . (9) 



2V2KC 
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The constant zq can be regarded as the position of the interface at t — 0, 
and Iq as its width. $ smoothly interpolates between the local minima of 
V: for z — > — oo and $ + for z — > +00. The corresponding values of the 
potential are 

^->=°- 

At $ m = SB/AC — $ + the potential V has a local maximum if A > 0. 
The substitutions Z — > — Z and t>o — > — t>o in formulas (7), (8) give another 
solution of Eq.(4), called anti-interface. 
It is clear that solution (7) exists if 

9B 2 > 32AC. (10) 

The parameter A has the following dependence on the temperature T [10|] 

A = a(T — T*), 

where a > 0. The constants 5, C and a, approximately do not depend on 
temperature. Condition (10) is satisfied if the temperature T is from the 
interval (T*, T c ), where T c is determined from the equation 9B 2 = 32aC(T c — 
T*). It is clear that T c > T*. For temperatures in this interval one phase is 
stable and the other one is metastable. 

The potential (2) can also lead to a static, symmetric domain wall. 
Namely, for the temperature T such that B 2 = 4AC the velocity v vanishes, 
V($ + ) = V($_) = 0, and the potential can be written in the following form 



V = C 

where for T = Tq 



($ - $ m ) 2 - $; 



fill 



In this particular case there is the degenerate ground state given by $ = $±. 
The potential (11) possesses the Z 2 symmetry 

$ -> 2$ m - $, 

and the interface becomes a static homogeneous, symmetric domain wall with 
the Z 2 topological charge. 



2.2 The comoving coordinates 

Here we quote the main definitions in order to introduce our notation. More 
detailed description of this change of coordinates, as well as a discussion of 
related mathematical questions, can be found in 0, ||], 0. 
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The comoving coordinates in the space, denoted by (a a ) = (a 1 , a 2 , a 3 = 
£) where a = 1,2, 3, are defined by the following formula 

x = X(a i ,t)+&(a i ,t). (12) 

Here x = where x a are the Cartesian coordinates in the space R 3 . 

The points X(a\t) form a smooth surface S which is parametrised by the 
two coordinates Q a 1 , a 2 . In general, the interface moves in space, hence X 
depends on the time t. 

The surface S is fastened to the interface - the shape of it mimics the 
shape of the interface and they move together. We shall see that for con- 
sistency of the perturbative scheme X(a\t) has to obey certain equation 
from which one can determine evolution of the surface S. The coordinate 
£ parametrises the axis perpendicular to the interface at the point X(a l ,t). 
The vector p(a\t) is a unit normal to S at this point, that is 

p 2 = l, pX tk (a\t) = 0, 

— * — * 

where X jk = dX /da . The surface S is characterized in particular by induced 
metric tensor on S 

9ik = X y iX^, 

and the extrinsic curvature coefficients 

K i:j = //.V , r 

The matrix (g tk ) is by definition the inverse of the matrix i.e. g lk g k i = 

si- 

The two by two matrix (Kik) is symmetric. Two eigenvalues fei, k 2 of the 
matrix where K* = g ll Kij, are called extrinsic curvatures of S at the 

point X. The main curvature radia are defined as Ri = 1/ki. Thus, by the 
definition 

Ki = ^ + ^, det(Ki) 



R\ R2 3 R1R2 

In general the curvature radia vary along S and with time. 

The coordinates (a a ) replace the Cartesian coordinates (x a ) in a vicinity 
of the interface. Components of metric tensor in the space transformed to 
the new coordinates are denoted by G Q/ 3. They are given by the following 
formulas 



G33 — 1, G^k — Gk3 — 0, Gik — N^irNj 



fcj 



x The Greek indices a,/3, ... have values 1,2,3 and they refer to the threedimensional 
space, while the Latin indices k, I, ... have values 1,2 and they refer to the inner coor- 
dinates a 1 , a 2 on the surface S. 
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where 

Dependence of G a p on the transverse coordinate £ is explicit, and a 1 , a 2 enter 
through the tensors g^ K\ which characterise geometry of the surface S. 
Components G a/3 of the inverse metric tensor have the form 

G 33 = x G 3k = G k3 = Q) Qik = ( N -l^ g rl( N -l}k^ 

where 

N 

and 



(N-'Yr = I [(i - ZKbsi + tK 



N = det(Nt) = 1 - iK\ + \i\K\K\ - K\K\) = (1 - - j-). 

In order to transform Eq.(3) to the comoving coordinates we use the 
standard formula 

A* = -LAfvgG^||V (13) 

where G = det{G a p), VG = y/gN, g = det(gik)- 

The time derivative in Eq.(3) is taken under the condition that all x a are 
constant. It is convenient to use time derivative taken at constant a a . They 
are related by the formula 

d d da 13 , d 

di lxa = di l * a + ^r lxa d^- (14) 

Finally, let us introduce the dimensionless variables s and (f) instead of, 
respectively, £ and $: 



£ = 2/ s, <S>£,ai,t) = -^JL < f>(s,<T i ,t). (15) 

The coordinate s gives the distance from S in the unit 2/ related to the 
width of the planar interface. 

Using formulas (13-15) we can write Eq.(3) in the following form, which 
is convenient for construction of the expansion in width: 

= IS + f + 2^o 2 ^v^ (^ViW,*) (16) 
-a0 + (1 + «)0 2 - 3 , 
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where 

is dimensionless transverse velocity of the surface S, the dot denotes the 
derivative d/dt\ a c, and 

K 

Formula (9) and the condition (10) imply that < a < 1 for temperatures 
in the range [T* , T c ] . 

The homogeneous planar interface (7) can be obtained from the evolution 
equation written in the form (16) in the following manner. As the surface 
S we take a plane, hence Kj = 0. Moreover, S is assumed to move with 
constant velocity v , hence 

pX — v — const. 

Finally, 

d<p \ -0 

because we look at the interface from the comoving reference frame, and 
(90/(9(7* = because of the homogeneity. Then equation (16) is reduced to 

<90 1 d 2 d> , 9 , , . 

The solution previously given by formulas (6-9) now has the form 

= 0o(s), v = a-^, (18) 

where 

exp(g - s ) 
1 + exp(s - s ) 

0o ( s ) smoothly interpolates between and 1. This corresponds to interpo- 
lation between the minima of the potential V if p is directed from 
negative towards positive z's. If we choose the opposite direction for p we 
obtain the anti-interface. The constant So corresponds to Zq from formula 
(6). 
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3 The expansion in width for curved interfaces 



Let us begin from a brief description of ideas underlying the calculations 
presented below. The set of solutions of the nonlinear, partial differential 
equation (3) is very large. We are interested here only in a rather special 
subset of it, consisting of solutions which represent evolution of a smooth 
interface. Moreover, even within this subclass we concentrate on rather spe- 
cial interfaces, called by us the 'basic' ones. Their defining feature is that 
one can find a comoving coordinate system in which the order parameter of 
the interface is essentially given by <j>o(s), formula (19), modified by small 
corrections which take into account the nonvanishing curvature. 

By writing the evolution equation in the form (16) we have shown that Iq 
can be regarded as a parameter analogous to a coupling constant - it appears 
in Eq.(16) only as a coefficient in several (but not all) terms. Therefore, one 
may hope that a systematic perturbative expansion in l will turn out useful, 
as it is the case with other perturbative expansions so numerous in theoretical 
physics. The perturbative series can be constructed in the standard manner: 
the sought for solution and the velocity v are written in the form 

= 00 + h<Pl + ^002 + -, V = Vo + IoVi + IqV2 + ---, (20) 

and inserted in Eq.(16). Coefficients in front of successive powers of l in this 
equation are equated to zero. Notice that after the rescaling £ = 2l s the 
expansion parameter l is present also in N and (iV -1 )^. In the zeroth order 
we obtain the following equation 

-^ = ^-«0o + (l + a)0^-4 (21) 

which formally coincides with Eq.(17). Therefore, we can immediately write 
the relevant solution 

1 , , is exp(s — C (a % , t)) , s 
v = a--, Ms,er\t)= y \ n/ + xv 22 

2 1 + exp(s - C (a l ,t)) 

There are however two differences between the planar solution (18), (19) 
and the solution (22). First, we do not assume homogeneity of the inter- 
face, therefore the constant s from formula (19) is replaced by the function 
C (a l ,t) of the indicated variables. Second, the surface S is not fixed yet, 
while in the former case it was a plane. 

It is convenient to rewrite Eq.(16) as equation for the corrections S<f>,Sv, 
which are defined by the formulas 

= o (s,a\t) + 5(f), v = v + Sv. (23) 
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After taking onto account the fact that </>o obeys Eq.(21) we obtain equation 
of the form 

L5<f> = /, (24) 

with the linear operator L 

i Id 2 , 1J 



+ (a a + 2 (a + 1 



l2 



1 dC 


1 9<t> \ 


l°" Q at 





and 



- 2J p(N-%g kr X r (X + 2l s'p) (5^ - C 0ji ^) 
- fef + (G^N (5<p, k - C , k °t)) 

+(30 o -a-l)(50) 2 + (50) 3 , (25) 

where vq and 4>o are given by formulas (22). Now it is easy to see that order 
by order in l we obtain enumerated by n = 1,2,... inhomogeneous linear 
differential equations of the form 

L(fi n = fn- (26) 

For example, 

A = (*J-*)l£- < 27 > 

In the whole perturbative scheme Eq.(21) is the only nonlinear equation for 
the contributions to the order parameter <fi. We show in the Appendix that 
Eqs.(26) can easily be solved with the help of standard methods - one can 
construct the relevant Green's function for L. It is remarkable that the same 
operator L appears in all equations (26), and that the form of it does not 
depend on the surface S. For these reasons calculation of the corrections 
(p n is reduced to relatively simple task of finding f n and calculating the one 
dimensional integrals over s shown in the Appendix. 

The perturbative Ansatz (20) and Eqs.(26) are two parts of the Hilbert- 
Chapman-Enskog method. The third and most crucial part consists of in- 
tegrability conditions for Eqs.(26) @. Such conditions appear because the 
operator L\ the Hermitean conjugate of L, has a normalizable eigenstate 
with the eigenvalue equal to zero. Such eigenstate is called the zero mode. 
Let us first find the zero mode for the operator L: inserting 0o m Eq.(21) 
and differentiating this equation with respect to s gives the following identity 

Lij r = 0, (28) 
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where 

as [1 + exp(s — G )J 
Notice that ip r exponentially vanishes for s — > ±oo. Because the opera- 
tor d/ds is anti-Hermitean with respect to the scalar product < gi\g% >= 
f*™ dsglg2, the operator L is not Hermitean. Its Hermitean conjugate has 
the form 

tt = ld#~( a ~bte~ a + 2(C " + 1)00 ~ 30 °' 
The operator L* has a zero mode too, namely 

LHi = 0, (30) 

where | 

^ = exp[(2a-l)(s-C )]Vr. (31) 

The function ipi vanishes exponentially for s — > ±oo because < a < 1 for 
all temperatures in the range (T*, T c ). For a = 1/2, that is when the interface 
becomes the domain wall, the two zero modes coincide. 

Let us multiply both sides of Eqs.(26) by ipi( s ) an d take the integral 
J-™ ds. The l.h.s. of the resulting formula vanishes because of (30), hence 

+oo 

ds inf n = (32) 

-00 

for n = 1,2, ... . It turns out that these conditions are nontrivial. In partic- 
ular, they give evolution equation for the surface S. 

It should be noted that the conditions (32) are in fact approximate, but 
the neglected terms are exponentially small. The point is that in order to 
obtain Eqs.(26) we use the expansions of the type 




1 - Uos/Ri 

(i = 1, 2), which are convergent for s < % where sm = min(i?i/2/ , R2/2I0). 
Therefore, when deriving conditions (32), the integration range should be 
restricted to \s\ < sm- Because of the exponential decrease of ipi and f n at 
large |s|, this will give exponentially small corrections to these conditions. 
We assume that the ratios R4/I0 are so large that we may neglect those 
corrections. 



2 The subscripts 1 and r stand for left and right, respectively. The point is that (30) can 
be written as ipiL = 0. 
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Let us compute the basic interface up to the order Iq. For n — 1, condition 
(32) gives 

Vl = -R + -R { ^ 



because 



a Q {a) 



/+oo 
dsipiip r = £(2a + l,3-2a) 
-oo 



does not vanish for a in the interval (0, 1). Here B denotes the Euler beta 
function. The function a (a) has the symmetric "U" shape in the interval 
a G [0, 1], with the minimum equal to 1/6 at a = 1/2, and the upper ends 
reaching the value 1/3 for a = and 1. The condition (33) implies that the 
surface S obeys the following evolution equation 

<y u 2a - 1 



It formally coincides with the well-known Allen-Cahn equation |TT 
Now Eq.(26) with n — 1 is reduced to 

L0i = 0. (35) 

It has the following solution which vanishes at s — >• ±oo 

^d^W, (36) 

where C\ is a smooth function of the indicated variables. 

For n = 2, we obtain from formula (25), after taking into account the 
results (33) and (36), 

f 2 = (2sK]K\ - v 2 )d s <p + (30o - a - l)Cf(d s <p ) 2 

+f + cf k X k kc Q ^ d s 4> + 2^=£, (g^y/gCwddo) . (37) 

Straightforward integration over s as in (32) can be a little bit cumbersome. 
This calculation can be significantly simplified with the help of the following 
identity 

2 / ds^[34> - (a + l)]d s( p Q( p n = [ dsdM n , (38) 



where in the case at hand n = 1. Identity (38) is obtained from Eq.(26) 
by differentiating both sides of it with respect to s, multiplying by ipi and 
integrating over s, just as in the derivation of the integrability conditions 
(32). The integrability condition gives 

a x {a)K)K{ -v 2 + % (-3A + g lk X, k XC , 
+2A 2 C + (2a - l)g lk C , t C 0tk + 2C^K)K\ = 0, (39) 
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where 



1 d ( _ d 



yfg da 1 \ da k 
is the Laplacian on the surface S and 

_ 1 dao(a) 



ai(a) = ao(a)~ 



da 



For a from the interval [0,1] the function di(a) is almost linear. In particular, 
Q (0) = -3,a(l/2) = 0,a(l) = 3. 

The integrability condition (39) leaves the freedom of choosing whether 
we keep nonvanishing Cq or v 2 . It is clear that we can not put to zero both 
of them unless KjKf = (then S is a plane). The choice Co = gives 

v 2 = a x {a)K)K\. (40) 



This implies a correction to the Allen-Cahn equation of the form 

^\{kp) = l 2 v 2 , 



K 

where 5i(.) denotes the first order correction to Xp. In consequence, there 

— * 

will be a correction to the solution X of the Allen-Cahn equation, and correc- 
tions to gik and K ik . These corrections have to be taken into account when 
calculating f k with k > 3. It is clear that this version of the perturbative 
scheme is rather cumbersome. 

On the other hand, if we put v 2 = 0, then the evolution of the surface 
S is still governed by the relatively simple Allen-Cahn equation (34). The 
integrability condition (39) is now saturated by the function Co - it has the 
form of evolution equation for Co, namely 



2L 
K 



d t Co - g lk d ak XXC ,i 



- A 2 C (41) 



-(a - \)g ik Co,iCo, k - C Q K)Kl = \a x {a)K\K\. 

Let us also check the third integrability condition. Formulas (20), (25) 
give 

+%g ir X r X(C 1 4> r ),i + %C 0>i (g ir X r pK ir X r X)si> r 
+2sK)K\d s ^ - 2A 2 (dVv) + AsKlj=(^ k C , k A) j 
-^(y/g^KlCoMj + ^(^gK^CoMj 

+2(30o -a- 1)0102 + 0?, (42) 
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where we have put V2 = 0. The term with <p2 contains a new function C2 
(see the Appendix), however this function will not appear in the integrability 
condition because the integration over s eliminates the term proportional to 
C 2 . This can be seen from the identity (38): on the r.h.s. of it we have f 2 
in which C2 is not present. We see already from formula (42) that we can 
choose whether to keep nonvanishing U3 or C±. For the same reason as in 
the case of n = 2 we choose U3 = 0, and the integrability is saturated by C\. 
The resulting evolution equation for C\ has the following form 



2L 
K 



dtd-^d^XXd^d 



A 2 d - K\K{C X 



-(a - ^dkCodjC, = - (fa + ai C + C 2 ) K^K)K{ - {K\f} 



+2 



gVdjKl-^idrXfl + K^idrXX)} diC (C + fa) (43) 



where 



d j (K* y /gd k C )(C + fa) - 2Ki k d k C Q djC Q \l + (2a - 1)(C + fa) 



a -< 

&2 — • 

a 



It is easy to proceed to the second and higher orders. Using formulas 
from the Appendix one can write general solution 0„ of Eq.(26). It contains 
the function C n ((j\t) which obeys an equation analogous to (41) or (43). 
This equation follows from the integrability condition (32) with n replaced 
by n + 2 if we put v n = 0. Due to identity (38), in the derivation of that 
equation we do not need explicit form of <f> n +i. In the present paper we will 
stop our considerations at the first order. 

In order to obtain a concrete basic interface solution we have to specify 
initial data for equations (34), (41), (43). There is no restriction on the initial 
data, except the obvious requirement that perturbative corrections of given 
order should be small in comparison with the ones from preceding orders. In 
particular, l Ci <^ 1,IqC 2 1 and Iq/Ri 1. If we consider the interface 
solution only up to the first order correction (36) then a simplification ap- 
pears: without any loss in generality we may adopt the homogeneous initial 
data 

C (a\t = t ) = 0, C 1 (a l ,t = t ) = 0. (44) 

This can be justified as follows. Local deformation of S by shifting a small 
piece of it along the direction p results in the corresponding shift of the 
coordinate s. Therefore for any given basic interface we can choose initial 
position of S such that Co(a\ t — to) — in formula (22) for <fto- This can be 
done at one instant of time, e.g., at the initial time. Values of Co and position 
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of S at later times are determined uniquely by Eqs.(34), (41), and in general 
C does not vanish. Notice however that such a shift wiil influence terms of 
the order /q and higher in formula (25) for / - due to the explicit presence 
of s in N, (iV _1 )k / is not invariant under the translations s — > s + Co- 

The r.h.s. of Eq.(41) vanishes for a — 1/2, that is in the domain wall 
case. In this case the initial condition (44) implies that Co = for all times, 
and in consequence Co disappears from the first order perturbative solution. 

As for Ci, the reason for the homogeneous initial condition is that the 
interface with the first order correction, that is 

= 0o(s-C o ) + /oCi^ r (s-Co), (45) 

can be regarded as 0o( s — Co + foCi) to the first order in l , so again we 
can cancel C\ at the initial time to by suitably correcting the initial position 
of the surface S. Let us stress again that this works only at the fixed time 
instant. 

The initial data (44) imply that at the initial time the order parameter 
is equal to 0o( s )- Hence, the only freedom in choosing initial data for 
we still have is position of the surface S. When the second and higher order 
corrections are included one has to allow for more general than (44) initial 
data for C and C±, nevertheless the initial form of always is uniquely fixed 
by these data and the initial position of the surface S. 

To recapitulate, the perturbative solution to the first order has the form 
(45) where 0o, are given by formulas (22), (29), respectively. Formula 
(45) gives dependence on s explicitly. The functions C , C\ are to be deter- 
mined from Eqs.(41),(43) with the initial data (44). One also has to solve 
Allen-Cahn equation (34). In certain cases these equations can be solved 
analytically, e. g., for spherical interface discussed in next Section. In gen- 
eral case one will be forced to use numerical methods. In comparison with 
the initial evolution equation (3) the gain is that the equations for S, Co, C\ 
involve only two spatial variables a 1 , a 2 . Such reduction of the number of 
independent variables is a great simplification for numerical calculations. 

The perturbative solution obtained above can be used in calculations 
of physical characteristics of interfaces. In the following Section we obtain 
formulas for local velocity and surface tension of the interface. In the Section 
5 we discuss evolution of a spherical interface. 
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4 Local velocity and surface tension of the in- 
terface 



Let us apply the expansion in width in order to find local transverse velocity 
and surface tension of the interface. We shall use the first order solution 
(45). The velocity is obtained from the condition <f> = const. It does not 

necessarily coincide with pX given by Allen-Cahn equation (34). Because we 
neglect terms of second and higher order in Iq, we may write <fi in the form 
4>o(s — C + l Ci) from which we see that is constant on surfaces given 
in the comoving coordinates by the condition s — Cq + IqC\ = s , where Sq 
is a constant. It follows from formula (12) that in the laboratory Cartesian 
coordinate frame these surfaces are given by x (a\t), where 

x (a\ t) = X(a\ t) + 2/ (s + Q, - kCi)p{a\ t). 

The transverse velocity of the interface is equal to x p. In order to calculate 
it we take time derivative of xo, project it on p, and use equations (34), (41), 
(43) with the initial data (44). The result can be written in the form 

7 ^ 2a - 1 1 1 , (I 1 \ , 2 / 1 1 \ tAn . 

The unit normal vector p is directed from the isotropic phase (<3> = $o) to 
the ordered phase ($ = $+). In the next Section we shall use formula (46) 
in the case of spherical droplets. 

The surface tension is another basic characteristics of the interface. It 
can be determined from a formula for free energy of the interface, which is 
defined as follows. The surface S cuts the total volume of the sample into 
two regions denoted below by I and II. Let us imagine that in region I there 
is the homogeneous isotropic phase with constant free energy density equal 
to V{$-) = 0, and in region II the homogeneous ordered phase for which 
V(& + ) = K^v /(961qC). The normal vector p points to region II. The free 
energy of the interface is defined as the difference 

Fi = F — V($ + )V n , 

where Vn denotes volume of the region II, and F is the total free energy of 
the sample given by formula (1). We shall compute Fj using the first order 
solution (45). Then, without any loss of generality we may put C = C\ = 
at the given time, as argued in the preceding Section, while the surface S 
remains arbitrary. Therefore, we need only 0o( s ) — exp(s)/(l + exp(s)). 
Because the dependence on the coordinate s = £/2/ is explicit, we can 
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integrate over s in the formula (1) for the free energy F. The volume element 
and the gradient free energy are taken in the form 

dx a dx a da a daP 

Neglecting terms quadratic in Iq/Ri we obtain the following formula for the 
free energy of the interface 



/ k dA, 
Js 



where 

K 

K 



l-£-2)(l-2a)(£ + £) 



(47) 



can be regarded as the local surface tension of the interface at the point 
a 1 , a 2 . -Ri,i?2 are the main curvature radia of the surface S at that point, 
and dA = ^fgdo x d<7 1 is the surface element of S. Of course, this formula for 
Fi can be trusted if lo/Ri <C 1. 

For a spherical droplet of the ordered phase embedded in the isotropic 
phase R11R2 are positive (the signs follow from formulas given in Section 
2.2) and of course equal to the radius of the sphere. If a < 1/2, formula 
(46) implies that the droplet grows (if its radius is large enough) because 
pxo < and p is the inward normal. In this case the curvature correction 
diminishes the surface tension, and k increases as the droplet grows. In the 
reverse situation - the isotropic phase inside and the ordered one outside - 
the curvature increases the surface tension and k decreases as the droplet 
grows. 

In the case of a growing droplet of radius R of the isotropic phase in the 
ordered medium k has the same dependence on the curvature. Here p is 
the outward normal, a > 1/2, px > 0, and R± — R 2 — —R. Nevertheless 
the values of surface tension in both cases are different because l present in 
formula (47) depends on a, namely Iq ~ (1 + a), see formula (57) below. 

Notice that the first order curvature correction to k vanishes in the domain 
wall case (a = 1/2). 



5 Evolution of spherical interface 

Let us now apply the formalism developed in Sections 3 and 4 to evolution of 
spherical droplets. We assume that a 7^ 1/2 in order to exclude the relatively 
simpler case of the domain wall. The surface S is parametrised by 

X = TR(t)p(9,ij). 
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Here 9, ifj are the spherical angles, and p is the inward (outward) normal to 

the sphere when a < 1/2 (a > 1/2). Thus, pX = =fR, with the upper sign 
corresponding to a < 1/2. The Allen-Cahn equation has the form 

iR k =k-h (48) 

where 

R 4/0 



|l-2a| 

Integration of equation (48) yields the following formula 

R(t) , ,R(t) 2K R(0) , . -R(O) 

— — + In | — — -1| = — = -£ + -^+ln| — — -1|. (49) 
R* R* -R*7 -R* -R* 

Evolution equation (41) for Co now has the from 

2^°° " RW Co = 2RW' m 
It has the following solution 

ai it>* fl(0) - iZ(t) 

Co(t) " i?(o)-^ (51) 

which obeys the initial condition C o (0) = 0. 

Evolution equation for Ci is obtained from the general equation (43) . For 
the spherical bubble it has the form 

T^Ci - j~ 2 Ci = T2(\a 2 + ai C + C 2 )-l (52) 

where as usual the upper sign corresponds to a < 1/2. We do not know 
explicit solution of this equation. 

There are two cases when the spherical droplets grow: a droplet of the 
ordered phase when a < 1/2, and a droplet of the disordered phase when 
a > 1/2. In both cases formula (46) for the radial velocity of expansion r 
gives 

7/0. 2a — 1 ^0.1 / \\A0\2 1 / C o\ 
^ r o = 4 R + ^ a ^W ~a 2 (a)(-y. (53) 

The expansion velocity is identical for all surfaces of constant 0. 

It is clear that there is a minimal R, let us denote by R m i n (&), such that 
ro > 0. We have found numerically that 

R 

R min( a ) = 3^T> ( 54 ) 
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where the function z(a) is symmetric with respect to a — 1/2 and it has 
values in the interval [0.866, 1.049]. For example, z(0) = 0.866, z(0.1) = 
0.984, z(0.20) = 1.045, z(0.3) = 1.040 z(0A) = 1.012, z(0.5) = 1.0. Notice 
that -R m j n (a) diverges when a — > 1/2. Thus, in the Ginzburg-Landau model 
nucleation of expanding droplets is possible only if we heat the ordered phase 
to a temperature above T , or cool the isotropic phase below T . 

For large time t, when the droplets are very large, the velocity r(t) be- 
comes equal to the velocity of the planar interface 

K 

r °c(a) = TT-rl 1 - 2a|, 
27/ 

as expected. Notice that -R mm («) and r^a) are not independent: 

2K 

z(a)f 00(a) R m [ n (a) = — . 

Parameter ct is related to the temperature: 

20 

a = 



1-29 + ^^19' 

where 



0= ~{T-%) 



8aC 
952 

is a reduced temperature. T and T c correspond to 9 = 2/9 and 1/4, cor- 
respondingly. The interval a G [0,1/2] corresponds to 9 e [0,2/9], and 
a G [1/2, 1] to G [2/9, 1/4]. The temperature dependence of l can be seen 
from formula 

,V2KC , x 

^o = (l + «)^-, (55) 

which follows from the definition (9) after some algebraic manipulations. 
i? mm is proportional to Iq/\1 — 2a\, which can be written in the form 



l V2KC 



\l-2a\ 6B \SVT^I9-1\' 

In the interval 9 G [0, 2/9] , which corresponds to T G [T*, 7b], Z /|l - 2a\ 
monotonically grows from V2KC/3B to infinity. 

Using formula (55) and the symmetry of z: z{\/2 — 5) = z{l/2 + 5), we 
obtain the following relation 



i? min (l/2 - j) = 3-25 
i? min (l/2 + 5) 3 + 25 



< 1, (56) 
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where 5 G (0, 1/2). Thus, the minimal size of the droplets of the isotropic 
phase which appear and grow when a > 1/2 is significantly larger than the 
size of the droplets of the ordered phase which can appear for a < 1/2. 

The velocities depend on temperature. In particular, comparing them 
for temperatures below and above T , 

7-^(1/2-5) = 3 + 25 i 
7^(1/2 + 5) 3-25 > 

Hence, the droplets of the isotropic phase expand more slowly than the 
droplets with = 1 inside. 

Our main goal in this paper has been to develop the perturbative expan- 
sion for the curved interfaces. We plan to apply it to interfaces in liquid 
crystals in a subsequent work. Nevertheless, just in order to get an idea 
what our formulas predict, we have estimated Iq and for interfaces in 
nematic liquid crystal MBBA. The model defined by formulas (1), (2) and 
Eq.(3) can be related to de Gennes-Landau theory in a single elastic con- 
stant approximation {L\ — K, L 2 — 0). We take data found in |12[ [13], [14]]: 
T* » 316K,a « 0.021J/(cm 3 • K),B w 0.07 J/cm 3 , C « 0.06J/cm 3 (after a 
change to our notation), and K ~ 6 ■ 10~ 12 iV, 7 ps* 5.2 ■ 10~ 2 kg/(m ■ s). We 
have identified 7 with the rescaled rotational viscosity jiLx/Ku at a temper- 
ature close to T*. Then, T — ps li^, T c — T* ps 1.2K. The width / and 
the velocity of the planar interface are given by the following formulas: 

o . . 1 1 — 2a I cm 

l w 40(1 + a) ■ 10 _8 cm, ^(a) w 1.4^ L — . 

1 + a s 

Notice that even for rather small droplets with the radius of several hundred 
Angstrom the ratio Iq/R is rather small. 



6 Remarks 

1. We have shown how one can systematically compute curvature corrections 
to the transverse profile and to the local velocity r of the interface. Due 
to the presence of functions Ct evolution equation for the surface S has 
the relatively simple form (34) and there are no curvature corrections to it. 
The formalism for interfaces is a generalisation of the one constructed for 
domain walls [^ |], |]. The main new ingredients are the C function and the 
(2a — 1)/Iq term in Allen-Cahn equation (34). By including them we have 
significantly enlarged the range of physical applications of the perturbative 
scheme. This justifies the present publication. 
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2. The model we have considered is special in the sense that the exact pla- 
nar interface solution 0o is known. Moreover, the solutions of the equations 
L<f)n = fn are given (almost) explicitly too, because the one dimensional in- 
tegrations in formula (60) below can be easily calculated numerically. Here 
the crucial point is that we know explicitly the two linearly independent so- 
lutions ip r and ip2 of the homogeneous equation Lip = 0. In other models, the 
analogs of 0o, VVj ^2 can be found at least numerically because the pertinent 
equations are relatively simple differential equations with the single indepen- 
dent variable £ (or s after a rescaling). In our perturbative scheme we need 
to perform explicitly only integrals over s, as in integrability conditions (32) 
or in formula (60) for <j> n . Such integrals can easily be calculated numerically 
also in the case when only numerical solutions </>o,Vv, ^2 are known. 

3. In our approach, in order to describe the evolution of curved interface 
we introduce the surface S, and the functions Ck which can be regarded as 
auxilliary fields defined on S and coupled to extrinsic curvatures of it. The 
corresponding evolution equations, that is Allen-Cahn equation (34)for S 
and equations (41), (43) and analogous ones for Ck, have one independent 
variable less than the original equation (3). This is significant simplification 
from the viewpoint of both computer simulations and analytical approaches. 
Therefore we think that our perturbative scheme is an interesting tool for 
investigations of dynamics of interfaces in Ginzburg-Landau effective models. 

4. The perturbative solution we have presented above is based on the planar 
homogeneous interface 4>o(s). Moreover, the dependence on the transverse 
coordinate s is uniquely fixed by the perturbative scheme once the initial 
position of the surface S and initial data for the functions Ck are fixed. This 
means that in our scheme we obtain a special class of interfaces, distinguished 
by the particular form of the dependence on s. In other words, the trans- 
verse profiles of the interfaces provided by the perturbative solution are not 
arbitrary. Intuitively, the interfaces can be regarded as the planar interface 
folded to a required shape at the initial instant of time, and modified by 
necessary curvature corrections. Therefore, it seems appropriate to regard 
the curved interfaces obtained in our paper as the basic ones. More general 
interfaces could be obtained by choosing more general initial data and solv- 
ing Eq.(3). For such generic interfaces no analytic perturbative approach is 
available. 

5. One of advantages of a systematic perturbative approach is that one can 
make trustworthy estimates of neglected contributions and in consequence to 
check whether a given perturbative result is reliable. Formalism with that 
level of control can be used in order to make straightforward predictions of 
dynamical behaviour of an interface, but perhaps more important application 
is to "inverse problems", that is determination of parameters of Ginzburg- 
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Landau effective theory. From dynamical behaviour of the interfaces one 
could reliably infer what values have parameters of the model. For example, 
for a given liquid crystal one could experimentally determine coefficients in 
front of higher powers of the order parameter like $ 5 , $ 6 or terms of the type 
§d a §d a § in the formula for free energy F. For a discussion of the form of 
F for liquid crystals see, e.g., [|T^] . 

One could try to generalize our perturbative scheme for calculating the 
curvature corrections to interfaces coupled to a noise. In that case the r.h.s. 
of Eq.(3) would contain an external stochastic force which in particular would 
lead to fluctuations of the planar interface. In our formalism the dependence 
of, e.g., the surface tension k on curvature radia Ri is explicit and it comes 
from purely geometric quantities like the metric tensor or Jacobian, while nu- 
merical coefficients in front of powers of lo/Ri are given by integrals over the 
s coordinate and they are determined essentially by properties of the planar 
interface. Therefore, one may expect that if the stochastic force is present, 
values of these integrals would have to be averaged over the stochastic en- 
semble. This is another interesting direction in which one could continue the 
present work. 



7 Appendix A. The equations L(j) n = f n 

In order to determine n we have to solve Eq.(26). Using standard methods 



151 it is not difficult to obtain appropriate solution. 
Let us shift the variable s, 

s = x + C , 

in order to remove function Cq from (f>o present in the operator L. Then 
Eq.(26) acquires the following form 

L4> n (x, a\ t) = f n (x, a\ t), (57) 

where 

Id 2 Id 
L = 2 dx 2 + ^ ~ 2^dx ~ a + 2 ( a + - 3<f>l(x), 

f n (x, a\ t) = f n (x + C , a\ t), 4> n (x, a 1 , t) = (p n (x + C , o\ t), 

and 

7 / n exp(x) 
<po{x) 



1 + exp(x) 



22 



In the first step we find two linearly independent solutions of the ho- 
mogeneous equation Lip = 0. The zero mode ip r (x) is one solution of this 
homogeneous equation, and the other one has the form 

l/j 2 ( X ) = 1p r (x)h(x), 

where ip r = d<fi /dx and 

h(x) = - j^+i exp[-(2a + l)x] - \2 exp(-2a:r) (58) 

+ i J k[ ex P[( 1 " 2a ) x \ ~ l \ + ih 2 ex P[ 2 (! - a ) x \ + 3^ ex P[( 3 " 2c 04 

For a — > 1/2 the first term in second line reduces to 6x. Having those 
solutions one can construct the relevant Green's function and the solution 

4>n, 

(f) n (x, a\ t) = -2ip r (x) Jq dyipi{y)h{y)f n {y, a\ t) 

+2M x )f-oodyMy)fn(y,v\t) + c n (a\t)i> r ( x ), (59) 

where 

Mv) = exp[(2a - l)y]^ r (y). 

The functions C n ,n = 1,2, are utilised to saturate the integrability con- 
ditions (32). 

The solution of Eq.(26) is given by the formula 

0„(s, a\ t) = 4> n (s - C , <t\ t). 

8 Appendix B. Stability of the interface 

Significance of our theoretical analysis of the interfaces depends on their 
stability with respect to small perturbations. It is sufficient to check the 
stability of the planar interface because our perturbative solution is based on 
it. Mathematically, the stability is related to the sign of eigenvalues of certain 
operator, and it is a model dependent property. Considerations presented 
below apply to the model defined by formulas (1-3), of course. In this case 
the linearised evolution equation for small amplitude perturbations S(f> of the 
planar interface 4>o(s), formula (19), in the comoving reference frame has the 
form 

^^0 = ^0 + 2/2(^ + ^)50, 

where L has been given below formula (24), and x l ,x 2 are the two Cartesian 
coordinates in the plane of the interface. Because L does not depend on 



23 



x 1 ^ 2 , we may pass to Fourier modes S(f> of S(f> with respect to these coordi- 
nates. After the substitution 

5(f) = exp[(- — a)s]\& 

we obtain the following equation 

K 

where 

N = -\d 2 s + \{a + l -f + 2l%k 2 - 2(a + 1)0 O + 30 2 , 

with k 2 = {k\) 2 + (fe) 2 , fci,A;2 being the wave numbers Fourier conjugate 
with x 1 , x 2 . 
Notice that 

= exp[(a - -)s]ip r , 

where ip r is the zero-mode given by formula (29), is eigenfunction of the 
Hermitean operator N. The corresponding eigenvalue is equal to zero. Be- 
cause \l/o does not vanish at any finite s, it represents "the ground state" of 
a fictitious system with iV as "the Hamiltonian". Hence, all other eigenval- 
ues of iV are strictly positive. Moreover, looking at "the potential" in "the 
Hamiltonian" N one can see that the zero eigenvalue and the next one are 
separated by a finite gap. The eigenmode corresponds to a parallel shift 
of the interface as a whole in the x 3 direction. All other eigenmodes decay 
exponentially with a characteristic time equal to t c = 21^ /(K A), where A de- 
notes the corresponding eigenvalue. Thus, the planar interface in our model 
is stable. Consequently, also the curved interfaces are stable with respect to 
small perturbations provided that their curvature radia are large enough for 
validity of our perturbative expansion. 

Let us point out that from point of view of applications in condensed 
matter physics even unstable interfaces can be interesting if unstable modes 
grow in time so slowly that the interface manages to travel across the sample 
before these modes become visible. The interface has finite normal velocity 
pxo, formula (46), and in any real experiment the sample occupies a finite 
volume. 
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